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We propose a scheme for dispersive readout of stored energy in one mode of a nonlinear super- 
conducting microwave ring resonator by detection of the frequency shift of a second mode coupled 
to the first via a Kerr nonlinearity. Symmetry is used to enhance the device responsivity while 
minimizing self nonlinearity of each mode. Assessment of the signal to noise ratio indicates that the 
scheme will function at the single photon level, allowing quantum non-demolition measurement of 
the photon number state of one mode. Experimental data from a simplified version of the device 
demonstrating the principle of operation are presented. 
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High quality factor superconducting microwave res- 
onators are versatile devices that promise to form an 
essential part of a scalable quantum computing architec- 
ture. They have been used to readout superconducting 
qubits (fltxx jlj, charge [2|, transmon Q), and to provide 
a long-lived memory for qubit states via "parking" 
Integrated with SQUIDs they have been used to imple- 
ment high-sensitivity parametric amplifiers Q , create 
tunable resonators for the selective coupling of multiple 
qubits 043) ^iid study deghasing due to Kerr coupling 
between resonator modes [10]. In two-qubit systems, res- 
onators function as a cavity bus to couple far separated 
qubits [ll| . Coupled to polar molecules [l3| and nanome- 
chanical resonators [l3l| they have been proposed for side- 
band cooling experiments and for the construction of ex- 
otic superpositions of states. Theoretical proposals have 
been advanced for construction of non-classical states by 
using the Kerr coupling between two nanomechanical res- 
onators and for dispersive readout of number states 
of these resonators using the anharmonic coupling be- 
tween two beam-bending modes [l^. Since supercon- 
ducting resonators may also be used as quantum mem- 
ory elements [16j], and as subcomponents of a qubit |i3i] or 
a quantum repeater (l7| . the ability to directly measure 
their quantum state is invaluable. 

In this paper we present a superconducting circuit that 
realizes a non-demolition measurement of the energy 
stored in a resonator mode. Realistic simulations and 
calculations based on measured noise data (isl . [l9| in- 
dicate that single-photon sensitivity will be achievable. 
Our proposal uses Kerr coupling between two modes of 
the resonator to measure the shift in the frequency of 
one mode as the other mode is populated. The non- 
demolition condition [l^l, that the state projected into 
by the measurement be an eigenstate of the system, is 
enforced in a precise way by parity symmetry in our de- 
vice, in which an odd-parity state is measured by an even- 
parity mode. Impressive progress on non-demolition 
measurements in superconducting devices, with single- 
quantum sensitivity, have been reported recently in var- 
ious qubit- resonator systems, with either the qubit [21] 
or the resonator photons |22h24| as the measured quanta. 
Our proposal is distinguished from this previous work in 




Figure 1: Proposed circuit for implementation of nondemoli- 
tion measurement using the Kerr effect. The superconducting 
ring resonator (red) is made from two transmission lines Ta,b 
of unequal lengths (as indicated) joined at the ends by non- 
linear elements realized by nominally identical DC SQUIDs 
Si, 2- The 180-degree hybrid causes the pump to couple only 
to ring modes with an odd voltage profile with respect to the 
horizontal midline of the device; the power splitter causes the 
probe to couple only to even modes. 



two ways: First, the non-demolition condition, which is 
generally violated by small off-resonant terms (as in [2^). 
is enforced exactly by symmetry in our circuit. Second, 
there is no qubit involved; both our measured and mea- 
suring system are weakly anharmonic resonator modes. 
A measurement scheme exploiting the Kerr mechanism 
between two such modes has also been invoked in [2^, 
but only as part of a scheme for measuring itinerant 
photons. The parameters necessary for achieving single- 
photon sensitivity in our circuit are readily accessible us- 
ing current superconducting technology, unlike the situ- 
ation in traditional quantum optics [2g. 

Figure [1] shows a schematic for the proposed circuit. The 
ring resonator circumference is X = 16.00 mm. The 
DC SQUID two-terminal response is that of a nonlinear 
inductance, described in leading order by 



Lsq — Lq 



(1) 



We use Lq = 0.19 nH and Li = 32 H/A^ f^, corre- 
sponding to an unbiased DC SQUID with critical cur- 
rent Id — 0.85 /iA per junction; inductance parameters 
can be tuned with flux biasing. For the transmission line 
segments T(j,b, we assume an inductance and capacitance 
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per unit length of £ = 0.45 pH//zm, C = 0.16 fF //^m, ap- 
propriate for coplanar waveguide (CPW) geometry with 
3 /xm center strip and 2 jiva. slots. The loss tangent and 
the effective dielectric constant are tan (5 = 10~^ and 
Cpt f = 6.45, typical for a high purity silicon substrate 
|27| . The ring resonator is excited using two frequency 
synthesizers (pump and probe) through four capacitors 
with C = 10 fF. 

Eq. |T]) is modified if we take account of the effect of 
another flux, that threading the large ring resonator 
loop, f&ring- This cffcct is vcry small because the self- 
inductance of the full ring Lt is quite large, ^ring slightly 
changes the DC operating point for the superconducting 
phase around the ring, thereby changing the phase values 
around which the AC dynamics occurs. The full expres- 
sion for the effective inductance including the effect of 

^rir^.n IS. 



27r V2Lo 



zLll 



(2) 



Here z = 27r<i>ri„g/$o. We estimate Lt = 16mm x 
0.45pH//im=7.2nH. For this value of Lt the modification 



of the leading term, and of the P term, is negligible. A 
new term proportional to / is present for z ^ 0, but its 
contribution will vanish in the rotating wave approxima- 
tion used below. 

The fundamental even and odd resonant modes, degen- 
erate in the ideal ring, are split in frequency by coupling 
to the SQUID inductances. In the linear regime the even 
(e) and odd (o) mode frequencies are 



/e,o = / - /./(I ± cos(7rQ!)) 



(3) 



where / 



and /,/ = pLoiZo with Zq = y/CjC. 



It will be experimentally convenient to explore cross-Kerr 
effects in two modes that are close in frequency; our odd 
and even modes can be split by any desired amount (from 
to hundreds of MHz) by choosing a, determined by the 
relative lengths of Tq and as shown in Fig. 1. We 
write a classical Hamiltonian for these two modes. The 
harmonic part is exactly as in while the additional 
anharmonic parts of the energy go like |Li(/e ± /q)'* (± 
for the SQUIDs 5*1 and 5*2 respectively), where /e/o are 
the mode currents at the location of the SQUIDs. 

In terms of the standard conjugate variables n and (j) Q 
we obtain a Hamiltonian describing two harmonic modes, 
with additional purely quartic nonlinear terms: 



1 C 2 r 

H = Ecenl + ELe^l + Econl + Elo4>1 + [Eie^t cos^ P + Eljl sin^ /3 -f 6El,El,4>l4'l sin^ /3 cos^ /?] (4) 



where /? — 7ra/2, and the parameters £'(c,L)(o,e)j the electric and magnetic energies in the modes, are as in [8]. Our 
symmetrical geometry causes many possible additional nonlinear terms, e.g., ones proportional to 4'e4'o^ to be absent 
[l8|). Eq. (HI can be quantized in the standard way in the rotating wave approximation (RWA) we obtain 

H^^ = hoJealae + huoalao + ^^^^^ 1^12^ El ^El^EceEc o [^ala^ala^ - aja^ - ala^] sin^ ^Scos^ /3-f 

El^Ece [{alaef - Qala,] cos^ p + EI^Eco [{o-Wf - 64aJ sin^ /^j. (5) 

I 



Here huji = 2^ EnEci- Note the cross-Kerr type term 
proportional to a\aeCLlao- Terms in the rotating frame 
that are proportional to aeOeaJaj^e''^^'^^"^'^")* (in the in- 
teraction picture) are dropped from the Hamiltonian. 
This term and other similar ones are not so rapidly time 
varying if /e ~ fo- For our circuit fe — fo < 50 MHz, 
which means that this additional term could definitely 
not be ignored in nsec-scale pulsed experiments [3]. How- 
ever, here we only consider the steady state circuit re- 
sponse. Further, we can, without changing any other 
aspect of the model, raise the even-odd splitting with 
simple stub-tuning techniques [l^. Thus, these terms 
need not be a concern in the discussion of the quantum 



aspects of our proposal. 

The circuit has been simulated using the Agilent Ad- 
vanced Design System(ADS) [1^ for values of a near 3/2, 
where the pump and probe nonlinearities are comparable 
(Eq. (U))). The resonance frequencies of the circuit are 
determined using linear analysis in which the circuit is ex- 
cited by sweeping the frequency of the pump and probe 
synthesizers and measuring the reflected response at each 
of the four ports. The reflected wave amplitudes are mea- 
sured as voltage drops Vi_4 across the output resistors 
R = 50 fl using the circulators to distinguish between the 
incident and reflected waves. The even and odd response 
is extracted as Mi_2 = Vi ±V2 (M3^4 = ¥^±¥4 could also 
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Figure 2: Simulated device response near a = 3/2, where q 
fixes the relative lengths of Ta and Th (see Fig. 1). Inset: 
phase response of M\ as pump power Ppump is swept from 
to 12 fW in steps of 2 fW for q = 3/2. fro,probc = 7.1322 GHz 
is the probe resonance frequency in the absence of any pump 
excitation, and Sfr.probe is the shift of the probe resonant 
response under probe excitation. 



be used). These low-power resonance frequencies match 
very well with Eq. ([3]). 

The finite-power circuit response is simulated using the 
Harmonic Balance (HB) tool in ADS [l^. This classical 
tool is suitable for extracting parameters for our simple 
quantum description, and for confirming that all signifi- 
cant aspects of the system's response is captured by our 
Hamiltonian. In these simulations, the pump synthesizer 
is set to the pump resonance frequency and the probe 
synthesizer frequency is swept to extract the probe reso- 
nance feature for different pump powers. The results of 
the HB simulation, for circuits with different values of a 
near 3/2, are shown in Fig. [S] 

From this data we extract the fractional shift of the probe 
frequency as a function of pump power. A cross Kerr 
nonlinearity is clearly seen, as expected from Eqs. (|4I5I) : 
at small pump powers, the frequency of the resonant re- 
sponse seen by the probe shifts linearly with pump power. 
The probe power is set to Pprobe — 5.3 fW so that we 
obtain a signal to noise ratio sufficient to detect single 
photons. At this probe power level the probe resonant 
response shows significant asymmetry due to the self- 
nonlinearity of the mode (Eq. ([S])); but it is still below 
the power where bifurcation takes place. At the highest 
pump powers it is evident that the circuit exhibits higher- 
order nonlinearities not described by Eq. ([5]). But the 
model is very good at the very low pump powers that are 
relevant for describing the behavior of the measurement 
when only a few photons are in the pump mode. 

We now calculate the expected signal to noise ratio, 
showing that single-photon detection is indeed possi- 
ble. The pump power at resonance required to maintain 
one-photon circulating in the resonator in this mode is 

Pphoton — '^^^fr,pump/Q^:P'f^'>^P^ where fr,pump: Qr,pump 

are its resonance frequency and quality factor, respec- 
tively. The quality factor of the pump mode of our 



circuit is coupling limited to ~3,000 and it has a reso- 
nance frequency of 7.12 GHz. Then Pphoton ~ 0.07 fW. 
By performing HB simulations of the phase response at 
Ppump — and 0.07 fW, we read off that the expected 
single-photon phase signal is 6sig — 19°. 

This signal level is to be compared with the rms devia- 
tions of the detected phase due to several sources of noise. 
Noise due to amplification is fundamentally no less than 
one-half a photon of power per signal mode, but in prac- 
tice is determined by an amplifier temperature T^r. The 
rms magnitude of the variations in the detected phase 
due to this noise is (l^ 

amp. 4:kBTnAl' 

^0 = • (6) 

For state-of-the-art HEMT amplifiers. Tat w 4K 
The optimal bandwidth is set by the probe linewidth, 
and is roughly estimated as Av « fr.probe/Qr,probe- Here 
Qr. probe IS the probe quality factor, which is also cou- 
pling limited to 3000, and fr.probe ~ 7.1GHz; thus, 
Av « IMHz. This gives 

^amp. _ 2.7° (7) 

SNR = 0stg/cre « 7.1 (single photon). (8) 

We need to confirm that amplifier noise is indeed the de- 
termining limitation on the signal to noise ratio. Since 
phase noise due to low-frequency fluctuations in the res- 
onant frequency is a dominant source of noise in ex- 
periments performed at higher power levels, we exam- 
ine this mechanism in detail. To estimate the effects 
of this noise, we extrapolate results from [30|. The 
fractional frequency noise of a CPW resonator similar 
to the one we have assumed is measured there to be 
Ssfii^o, Po)/fr = lO-^VHz at bandwidth vq = 1 kHz. 
Here Pq = — 60 dBm = 1 nW corresponds to what we call 
probe power. Frequency noise power can be converted to 
phase noise using S'e(i^o, ^o) = ^Ql,probeSsfM/ fr = 
3.6 X 10"" rad^Hz. The phase noise in the bandwidth 
Ai' of our system can be estimated as Sg{Ai^, Pq) — 
{Av/v^y^'^ Se{vQ, Po) = 1.1 X 10-12 radVHz. The 
resonator phase noise is drive power dependent and needs 
to be scaled for the probe power used in our simu- 
lation. For a coupling-limited resonator, the internal 
power is given by P = PprobeQr, probe', then the phase 
noise power for our parameters can be estimated to be 

S0{v, P) = {P/Po)~^^'^ Se{v, Po) = 8.3x10-12 radVHz. 
Finally, the phase measurement error is: 

af- = ^Se{v, P)/t- 0.18° (9) 

Thus, frequency fluctuations make an insigniflcant con- 
tribution to the SNR. We are unaware of any other mech- 
anisms that are likely to contribute at this level. 

Figure [3] shows that our value of SNR=7.1 will be suf- 
ficient to discern the quantum nature of the resonator 
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state. We plot the probability density function for re- 
peated measurements of the Kerr-induced phase shift 
0. in the probe microwave signal. This is given by 
a sum over the possible number-state outcomes, con- 
volved with a function describing the additive noise in the 
measurement: P{0) = X^j^o ^(^l"-)^!")' where P{n) 
is the probability of projecting the coherent state into 
the n*'' number state given by the Poisson distribution 
(P(n) = for Uav = 1), and P{9\n) is the probabil- 
ity of measuring the given that the mode has been 
projected into photon number n. We assume Gaussian 

additive noise P(6\n) — — exp ( — I where 

ag and 9sig are estimated above. The x-axis of Fig. 3 
is marked with the phase shift nO^ig corresponding to 
the different photon numbers n in a noiseless measure- 
ment. We see that even for SNR~ 3, much less than the 
estimated value, our circuit should readily detect quanti- 
zation of the coherent state. Various nonidealities, such 
as non-identical SQUID critical currents (Jc), introduce 
extra nonlinear terms which reduce our response. How- 
ever, our simulation shows that expected ±5% variations 
of Ic reduce the response only slightly. In fact, there is 
considerable scope for improving the SNR by simultane- 
ously measuring amplitude and phase response and 
by carefully designing the resonator geometry to reduce 
phase noise [l9j . 

We have performed an experiment which shows that the 
basic effect predicted here is readily seen, even in a de- 
vice that is not carefully designed to optimize the cross- 
Kerr effect. The circuit is a linear niobium transmission 
line resonator with a quarter-wave geometry, capacitively 
coupled to a readout line at the top (Fig. ^ and shorted 
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Figure 3: Predicted histogram of phase measurement out- 
comes 9 (measured relative to the zero pump power phase) 
for different measurement signal to noise ratios (SNR) for the 
device in Fig. 1. Measurement is simulated for a coherent 
state with average photon number Uav = 1- The quantum 
nature of the radiation field is not evident for SNR=1, but 
for the predicted values of SNR=7.1 clearly resolved photon- 
number peaks will be seen. 
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Figure 4: Experimentally measured cross-Kerr response of 
simple SQUID-tuned resonator. The circuit schematic is 
much simpler than the optimized one above (Fig. 1). A CPW 
resonator with 3 /xm center strip and 2 pm slots was used. The 
DC SQUID loop had an area of 5x5 jsce? with each junction 
having Ic = 1.2 /xA. The SQUID was not flux biased for data 
in this figure. Readout line was used for pump excitation as 
well as measuring the forward scattering parameter 521 for 
the probe signal. Inset shows resonant phase response of the 
probe for different values of the pump power. Probe reso- 
nance frequencies are plotted in the main plot; the blue line 
is guide to the eye. The initial linear behavior of this curve 
at low pump power (slope indicated by black dashed line) is 
a manifestation of the cross-Kerr effect. 

to ground through a nonlinear inductor (an Al/A10a;/Al 
DC SQUID) at the bottom end. This device was mea- 
sured in a He-3 cryostat at 360 mK. 

It has fundamental (probe) and third harmonic (pump) 
resonances at 2.68 GHz and 8.05 GHz, respectively. The 
probe mode has a coupling-limited quality factor Qr — 
2000. The upper right inset of Fig. |4] shows the phase 
response of the device for different pump drives. From 
these we obtain Kerr shift of the probe as a function of 
the pump mode power. The resonance frequency is eval- 
uated as the frequency at which the derivative of phase 
response with respect to the readout frequency is max- 
imum. The probe power is held constant at -115 dBm. 
A detectable Kerr shift is observable for pump power as 
low as -70 dBm. A decrease in resonator quality factor 
similar to [To'| is observed as well. While our choice of ma- 
terials and operating temperatures allows us to measure 
these nonlinearity-induced effects at relatively low pump 
powers, accessing the quantum regime requires carefully 
optimised design proposed in Fig. [T] 

To summarize, we have presented and analyzed a cir- 
cuit for dispersive readout of energy stored in one mode 
of a superconducting microwave resonator by measuring 
the frequency shift of another mode coupled to the first 
through a cross-Kerr nonlinearity. The proposed device 
is sensitive enough to operate at the single photon en- 
ergy level, and it can be used to perform a quantum non- 
demolition measurement of the photon number state of 
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the resonator. The basic effect is readily seen in a simple 
experiment using an unoptimized version of the device. 
Many other uses can be envisioned for this device; for ex- 
ample, modulation of the SQUID flux bias at — fo will 
implement a beam splitter between the two modes [28J . 
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